Entanglement control in coupled two-mode boson systems 
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A two-mode boson model, widely used for the physics of fast rotating nuclei and Bose-Einstein 
condensates, is studied in the context of entanglement control. We derive an analytical expression 
for the entanglement between the fields in this model as a function of time. We found that depending 
on the interaction strengths between boson modes and the nature of the initial boson states the 
dynamical evolution of the entanglement and the squeezing can occur independently. 
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There is currently an enormous effort underway to 
understand the rich dynamics of interacting Bose sys- 
tems. Among various problems are the properties of 
Bose-Einstein condensate (BEC) in external fields pj, the 
Bose coherent effects of excitons and polaritons in semi- 
conductor microcavities 0] , and the production of scat- 
tered radiation due to interaction of the incident laser 
beam with vibrational modes of a medium (the Raman 
effect). The generic feature of interacting Bose systems 
is the formation of Bose-field collective states with non- 
classical (squeezed, sub-Poissonian, etc) statistical and 
fluctuation properties. Another prominent feature is an 
entanglement produced by the interaction between differ- 
ent system constituents, which is one of the most subtle 
and intriguing phenomena in nature (see for a recent re- 
view 0,I3)- Nowadays, there is explosive activity in the 
study of the entanglement due to its potential usefulness 
in quantum teleportation, quantum cryptography, and, 
in general, in quantum information theory. 

The dynamical interplay between quantum entangle- 
ment and nonclassical properties of various Bose systems 
can be traced within a model of two-coupled harmonic 
oscillators. A particular example of interest is the Hamil- 
tonian 



H 



+ihgi (c\c 2 



ihg 2 (c\( 
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Here Ci )2 (c\ 2 ) are the annihilation (creation) boson op- 
erators of the fields with the energies tvjj\ t2 . The coupling 
between these fields is governed by dint of the coupling 
constants 31,2- Note that the bilinear form of the Hamil- 
tonian ([T]) corresponds to a linearized version of some 
more general interactions. 

This model has been applied in nuclear physics [f| 
and for a rotating BEC (cf. Ref. 0). In condensed mat- 
ter physics the Hamiltonian ([TJ) is used: i) to study the 
interaction between an atom and a radiative field Q; 
ii) as a starting point for analysis of electronic proper- 
ties of two-dimensional quantum dots in a perpendicular 
magnetic field [§, |T(| • The dynamics of two-component 
Bose condensate trapped in a double-well potential can 
be also mapped on the time-evolution of two coupled- 
harmonic oscillators in the low excitation regime [ll| . In 



the simplest case, the model describes two levels of the 
condensed atoms which are coupled owing to the clas- 
sical field of radiation [12]. In this case the interaction 
constants (<?i and 172) in the Hamiltonian |T]) simulate 
the coupling between the electromagnetic field and the 
two-level (two-mode) boson system. Varying the inter- 
action (the coupling) and the choice of initial states one 
can analyze various properties of this system. 

The bilinear form of boson interaction enables one to 
investigate the system dynamics analytically. In par- 
ticular, the exact solutions as well as the generation of 
squeezed states in this model have been analyzed in [l3| . 
The main aim of the present paper is to gain a better in- 
sight into the dynamical interplay between the strength 
of the interaction, entanglement and squeezing dynam- 
ics of the system. Here, we focus upon the dynamics of 
Gaussian states, which are of great practical importance. 

By means of the Bogoliubov transformation 
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the Hamiltonian Jl} is reduced to the diagonal form H 
Y^t=i {o\flk + 1/2^ with the two eigenmodes 



nl = [ V+ + 2(g 2 1 -gl) + (-l) k+1 A}/2 : 
A = [rt + Arj+igl - g\) + 8^ 2 (.g 2 + g z 2 )] 



(3) 
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where i]± = w 2 ±w 2 . The coefficients A m and B m of the 
transformation matrices in Eq. are defined in Ref. [T3 . 

By virtue of the inverse transformation and time evolu- 
tion of the collective operators ak{t) = a/c(0)e 



—iflkt 



one 



can determine the time dependence of the initial opera- 
tors c m 

2 

Cm(t) = E hnn(t)c n (0) + /3 m „(i)4(0)] , (4) 



n=l 



where 



k=l 
2 

P mn (t) = ]T [A^B k e~ m ^ - B k m A* k ^] . (5) 
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The matrix elements a mn (t), j3 mn {t) obey the relation 



n=l 



[\a mn {t)\ 2 -\p mn {t)\ 2 } =1. 



(6) 



In the most general case it is natural to expect that at 
the initial stage of the time evolution the fields can be 
found in the superposition of coherent and chaotic states. 
For example, the Bose condensate of "cold" atoms can 
be described by a coherent state. However, there is al- 
ways a nonzero temperature which spoils the plain coher- 
ence and brings some chaoticity (the decoherence effects) . 
Similar phenomenon occurs at the interaction of the laser 
beam with vibrational modes in media, which are char- 
acterized by some thermal distribution. Therefore, we 
assume that the initial density matrix can be presented 
in the following factorized form 



m = n 



(l + (^}) fc 
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In Eq.([7]) we introduced new operators &i j2 — c\ 2 — ci!i,2) 
where ol\ %2 are initial coherent amplitudes of the fields ci j2 
and the averages (rij) are associated with mean numbers 
of the bosons in the corresponding chaotic states. 

To trace the time evolution of the quantum state of 
the fields c\ and c± governed by the Hamiltonian JTJ), it 
is convenient to introduce the Wigner function, instead 
of the density matrix. It can be done by dint of the 
symmetric characteristic function 



X(M1> M2! t) = Tr I g(0) exp (/i m c]„(i) - n* m c rn (t)) 

{ \_m=l 

(8) 

It is straightforward to show with the aid of Eq. (T4|) that 
the Wigner function obeys the relation 

W{v x ,v 2 ;t) = —g J d 2 /iid 2 M2X(Mi,M2;i)expI? 



= W(u 1 (t) > u 2 (t);Q). 



(9) 



Here V = Em=l( !/ m/ i m ~ "mfe)i v i = 
ELifcW^W - Vm{t)Vm(t% and £ m (t) 
ELi [<m(*)£« - /3nm(*)C] (where £ = n or v). It 
means that for initial Gaussian states the Wigner 
function retains its Gaussian distribution in the process 
of evolution being specified by the Hamiltonian |T|) . 

Thus, if the fields c\ and c± are initially in the states 
being superpositions of coherent and chaotic states, i.e., 
their Wigner functions have the form 



I 
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(10) 



the system Wigner function at time t can be described 
by the following expression 



W( Cl ,c 2 ;t) = exp{- [T! + T 2 }} /Tr 2 Af (II) 



with 



N = M x M 2 = «m) + 1/2) x ((na) + 1/2) , 
T x = A^t)^ 2 + A 2 {t)\c 2 \ 2 - {Bl{t)c\ 

+ B*{t)c 2 + B*(t)c lC2 - Bl{t)c lC * 2 + c.c.) , (12) 



T 2 = Y j {\a l \ 2 IM l + [Dt(t)a + c.c.]) . 



(13) 



Here, the coefficients A, B, D are defined as 

Ki(i)| 2 + |/3ii(i)| 2 , \a 12 (t)\ 2 + \p 12 (t)\ 2 



Mt) = 
Mt) = 

Bi(t) = 
B 2 (t) = 
Bs(t) = 
+ 
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The Gaussian character of the system state (fTTj) is well 
suited to introduce the en tang lement measure within the 
approach proposed in RefQjl- The measure of entangle- 
ment between the fields can be calculated analytically in 
the form of the logarithmic negativity via the symplectic 
spectrum of the partial transpose of the covariance ma- 
trix. To proceed along this line we present the bilinear in 
fields part of the Wigner function (fTTj) by dint of the re- 
lation T\ — \C, T 'V Here, the inverse variance matrix 
of the form 



j I Xi Y 
X 1 Y T X 2 



(15) 
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is determined by the matrices 

X,; = 



Ai - 2RcB l -2InxBj 
-21ms, Ai + 2RcB t 



Y = 



Re(B 4 - B 3 ) -lm(B 4 + B 3 ) 
-lm(B 3 - B 4 ) Re(B 4 + B 3 ) 



i= 1, 2, (16) 
(17) 



And £ T = (qi,pi,q2,P2) is a transposed four-vector 
whose elements are the quadrature-component variables 
defined as Cj = (qt + ipi)/y2. 




FIG. 1: Entanglement Eq. (|18|) between the modes being ini- 
tially in coherent states. Model parameters are o>i = 1, u>2 = 2 
and (a) gi = Q2 — 0.5, (b) <?i = 0, (?2 = 0.5. 



As a result, the logarithmic negativity 
1 



E 



log 2 (4A) 



where 



.4 = 6- V^ 2 -detV, 

_ det Xi + det X 2 - 
B = det Y , 



(18) 

(19) 
(20) 



determines the strength of the entanglement for E > 0. 
Here, Xi and Y are diagonal and nondiagonal block ma- 
trices of the variance matrix V, respectively. For E < 
the composite state is separable according to the Peres- 
Horodecki criterion [l5| (see also Refs.flo]). This mea- 
sure enables us to trace the evolution of non-classical 
correlations in the system and the conditions for their 
amplification. 

The considered model provides various possibilities to 
study cumulative effects produced by a different strength 
of coupling between the fields prepared in the initial 
states of different degree of chaoticity (different temper- 
atures). In the present paper we discuss a few interest- 
ing cases of the dynamical control of the entanglement 



by dint of different choices of the initial states and cou- 
pling constants. For illustration we consider the case 
u>i = 1, u)2 = 2 (in relative units). Note, that at a 
given energy fkui of the c^-field initial state, the tem- 
perature Ti fixes the mean number of bosons (rii), i.e., 
(m) = {exp [fUvi/(kQTi)] — 1}~ . All frequencies and 
coupling constants are determined in units of the fre- 
quency u)\ and the time is scaled by inversed u)\. 

At zero temperature, when both the fields are in co- 
herent (vacuum) initial states, the interaction with equal 
coupling constants (gi = g^) produces the entangled sys- 
tem state (see Fig. la). However, the degree of entan- 
glement is not monotonic with a time. It displays an os- 
cillatory character with the distinctive alteration of the 
entanglement maxima and minima that can be consid- 
ered to be revivals and collapses. In this regime the fields 
ci and £2 were shown to be evolved into squeezed states 
[l3| that is a consequence of such quantum correlations. 
In case the system is associated with the interaction of a 
light with ion oscillations in a crystal (optical phonons), 
a suppression of zero-point optical vibrations (squeezing) 
may affect the electron-pnonon interaction in nanostruc- 
tures. We speculate that the nonmonotonic behavior of 
the entanglement and the squeezing results in the inter- 
action strength which would fluctuate in time. Note that 
this interaction being essential to the understanding of 
electron-spin decoherence effects in quantum dots is as- 
sumed to be independent on time (cf [TtJ )• 
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FIG. 2: The same as in Fig. [T] for (a) chaotic states with 
(m) = (712) — 0.5 and (b) coherent and chaotic states when 
<n 2 > =5. 

The entanglement also occurs for the interaction 
regime gi — and 32 7^ (k^T = 0). In this case the 
interaction produces the picket-fence effect for the entan- 
glement of equal maxima without revivals (see Fig. lb) 
but with the periodical alteration of its minima. How- 
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ever, in this regime the fields Si and ci remain in coherent 
states, i.e., the squeezing does not take place |13| . It ap- 
pears that the interaction that creates (annihilates) the 
boson pair (associated with the strength $2) is enough 
to produce the entanglement. On the other hand, both 
interactions (associated with the strengths gi and 52) are 
indispensable in the Hamiltonian ([T]) in order to produce 
the squeezing of the coherent initial states. The inter- 
action of the gi 7^ and 52 = type does not yield 
an entangled state for the system at all, as well as the 
single-mode squeezing for the coherent states. 

Thermal fluctuations are expected to attenuate the en- 
tanglement. In particular, when both the fields are ini- 
tially in chaotic states, the degree of their entanglement 
in the regime of equal coupling constants becomes very 
small (see Fig. 2a). Moreover, there is a critical temper- 
ature above which the entanglement disappears and the 
system becomes separable (for equal coupling constants 
in our choice of the initial energies w lj2 it corresponds 
to (711,2) ~ 0.63). However, if one of the initial states is 
found in a coherent state while the other is in a chaotic 
one, the entanglement evolves at regular intervals with 



a periodical alteration of maxima (see Fig. 2b). With a 
proper combination of the initial energies of the system 
and a choice of number of bosons in the chaotic (initial) 
state, one may maintain the entanglement at a noticeable 
level, even at high temperatures. 

In conclusion, the considered two-mode boson model 
demonstrates that the entanglement can be controlled in 
time by appropriate choice of the interaction at different 
temperature regimes. The obtained results show that 
this plain system exhibits a rich dynamics from the view 
point of the quantum information theory and its possible 
physical applications. 
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